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CHAPTER  20 


The  Significance  of  Sediment-Flow  Dynamics  on  Clay  Microstructure 
Development:  Riverine  and  Continental  Shelf  Environments 


Huon  Li  and  Richard  H.  Bennett 


Introduction 

The  highly  variable  dynamics  of  coastal  and  continental  shelf 
waters  have  a  significant  effect  on  the  transport  and  deposition  of 
sedimentary  materials.  The  complex  sediment  transport  mech¬ 
anisms  in  continental  shelf  waters  involve  the  interactions  of 
multiphase  flows  of  solids,  liquid,  and  organics  (dissolved  and 
particulate)  with  multicomponent  chemical  and  minerals  spe¬ 
cies.  Also  biological  processes  are  pervasive.  Fine  particles  are 
transported  as  individual  particles  and  floccules  of  various  sizes 
and  they  interact  with  the  fluid  flow.  The  solid  particle  coagula¬ 
tion,  breakup,  deposition,  erosion,  and  transport  processes  and 
mechanisms  depend  on  the  tidal  cycles  as  well  as  the  waves  and 
currents.  Particle  associations  at  the  sediment  water  interface 
are  directly  related  to  the  particle  associations  in  the  suspension, 
and  thus  the  microstructure  at  the  sea  floor  is  considered  to  be 
related  to  the  floe  sizes  in  the  water  column.  The  particles  (indi¬ 
vidual  particles  and  floes)  in  suspension  that  are  deposited  even¬ 
tually  at  the  depositional  interface  are  the  fundamental  “building 
blocks”  of  the  bottom  sediment. 

The  objectives  of  this  chapter  are  to  address  (1)  the  signifi¬ 
cance  and  role  of  clay  microstructure  in  flow  dynamics  (empha¬ 
sis  is  on  clay  microstructure  interactions  with  tide,  wave,  cur¬ 
rent,  and  turbulence),  (2)  the  importance  of  scale  in  relation  to 
microstructure  and  flow  dynamics,  (3)  the  role  of  particle-to- 
particle  interaction,  bonding,  and  shear  that  occurs  near  the 
sediment-water  interface,  and  (4)  the  role  of  the  “outer  region” 
in  sediment  transport  and  microstructure  development.  Empha¬ 
sis  is  on  the  physicochemical  interactions  rather  than  biochemi¬ 
cal  or  biological  processes,  although  these  important  processes 
can  not  be  neglected  in  the  final  analysis  microfabric  develop¬ 
ment.  A  comprehensive  treatment  of  the  field  of  sediment  trans¬ 


port  and  dynamics  is  beyond  the  scope  of  this  chapter  and  thus 
limited  reference  material  is  discussed. 

The  surficial  sediments  are  often  reworked  by  organisms 
(bioturbation)  that  modify  the  physical  and  mechanical  properties 
and  microstructure  of  the  sediment  (Richardson  and  Young,  1980; 
Bennett  and  Nelsen,  1983).  The  clay  microstructure  of  continen¬ 
tal  shelf  sediment  has  been  studied  extensively  by  Bowles  ( 1 968a, 
b),  Bennett  (1976),  and  Bennett  et  al.  (1977,  1979,  1981).  Com¬ 
prehensive  studies  of  the  clay  microstructure  and  the  physico- 
chemistry  of  clays  have  been  presented  by  Grim  (1968),  Mitchell 
(1976),  and  Bennett  and  Hulbert  (1986).  Collective  papers  on 
shelf  sediment  transport  can  be  found  in  the  book  edited  by  Swift 
et  al.  (1972)  and  papers  on  cohesive  sediment  dynamics  appear  in 
the  lecture  notes  edited  by  Mehta  (1986). 


Clay-Water  Interactions 
Clay  Mineral  Particles 

Clay  minerals  are  crystalline  materials  that  are  uniquely  layered 
in  plate-like  structures.  Clay  minerals  are  classified  as  hydrous 
aluminum  silicates  (phyllosilicates).  Their  structures  are  com¬ 
monly  composed  of  composite  layers  built  of  components  of 
tetrahedrally  and  octahedrally  coordinated  polyhedra.  The  very 
common  clay  minerals  are  kaolinite,  smectite,  illite,  and  mixed- 
layer  clays.  The  clay  mineral  particles  are  usually  less  than  2  pm. 
Details  of  clay  mineralogy  and  classification  can  be  found  in 
works  by  Grim  (1968)  and  Mitchell  (1976).  Selected  charac¬ 
teristics  of  the  common  clay  mineral  particles  are  listed  in  Table 
20. 1 .  The  shape,  size,  cation  exchange  capacity,  specific  gravity, 
and  specific  surface  are  important  properties  of  the  clay  minerals 
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Ibble  20.1.  Characteristics  of  clay  mineral  particles.* 


Mineral 

Shape 

Size 

(um) 

Cation 
exchange 
capacity 
(meq/100  g) 

Specific 

gravity 

Specific 

surface 

(m2/g) 

Kaolinite 

Six-sided 

flakes 

0.1-4  x 
0.05-2 

3-15 

2.60-2.68 

10-20 

Montmoril¬ 

lonite 

(Smectite) 

Flakes  (equi- 
dimensional) 

>  10A  x 
up  to  10 

80-150 

2.35-2.7 

50-120 

primary 

700-840 

secondary 

Illites 

0.003-0. 1 
X  up 
to  10 

10-40 

2. 6-3.0 

65-100 

♦From  Mitchell  (1976). 


that  are  important  in  clay  microstructure-flow  interaction 
behavior.  For  example,  floculation  and  settling  velocity  are 
affected  by  the  shape,  size,  specific  gravity,  cation  exchange 
capacity,  and  specific  surface  of  the  interacting  particles.  Atoms 
of  clay  minerals  are  bonded  into  a  three-dimensional  structure, 
and  the  termination  of  this  structure  at  a  surface  produces 
unbalanced  force  fields,  positive  or  negative  charges  depending 
on  the  nature  of  the  termination  or  break  (Bennett  and  Hulbert, 
1986).  Isomorphous  substitution  within  the  crystal  lattice  also 
gives  the  clay  particle  a  net  negative  charge.  When  clay  mineral 
particles  are  2  pm  or  less  in  size,  the  surface  force  fields  become 
a  significant  factor  in  clay  microstructure  development.  Due  to 
the  platy  shape  of  the  clay  mineral  particles  and  the  fact  that  clay 
surfaces  in  general  have  a  residual  negative  electrostatic  charge, 
surface  and  colloidal  forces  exert  important  influences  on 
particle-to-particle  interactions.  To  preserve  electrical  neutral¬ 
ity,  exchangeable  cations  are  attracted  and  held  on  the  surfaces 
and  the  edges,  and  in  some  clays  between  the  unit  cells.  The 
quantity  of  these  exchangeable  cations  necessary  to  balance  the 
charge  deficiency  of  a  particular  clay  mineral  is  called  the  cation 
exchange  capacity  (CEC).  The  CEC  is  expressed  as  milliequiva- 
lents  per  100  g  of  dry  clay  (Grim,  1968;  Mitchell,  1976). 


Suspended  Clay  Particles— Salt  Water  System 

The  behavior  of  the  suspended  clay  particles  in  salt  water  is 
extremely  complex.  The  flocculation  behavior  of  clays  in  sus¬ 
pension  is  dependent  on  several  chemical,  physical,  and 
environmental  properties  such  as  electrolyte  concentration,  die¬ 
lectric  constant,  temperature,  ion  size,  adsorbed  water,  pH,  and 
clay  mineralogy  (Bennett  and  Hulbert,  1986).  The  particles  may 
coagulate  or  remain  stable  depending  on  the  balance  between  the 
electrostatic  forces  and  the  van  der  Waals  forces  and  environ¬ 
mental  conditions.  The  presence  of  organic  matter  is  an  impor¬ 
tant  factor  in  the  aggregation  of  clay  particles.  In  addition,  other 
complex  variables  include  the  bonding  of  clay  particles  by  mucus 


Tible  20.2.  Characteristics  of  two  mud  groups  according  to  Ohtsubo  and 
Muraoka  (1986). 


First  mud  group 

Second  mud  group 

Type  of 

Kaolinite 

Bentonite 

cohesive 

Kaolinite  +  organic  matter 

Bentonite  +  organic  matter 

sediment 

Montmorillonite 

Natural  muds 

Exchangeable 

cation 

A3*,  Ca2*,  H* 

Na*,  Li* 

Resuspension 

1 .  Clear  interface  between 

1 .  Obscure  interface 

behavior 

bed  surface  and  upper 

between  bed  surface  and 

water 

upper  water 

2.  Cracks  in  the  bed  surface 

2.  Hazy  streamwise  streak- 

just  before  bed  destruc- 

lines  consisting  of  tiny 

lion 

sediments 

3.  Bed  destruction  begins  at 

3.  No  cracks  in  the  bed  sur- 

the  cracks 

face  for  any  current  con¬ 
dition 

4.  No  flow  of  bed  at  bed 

4.  Flow  of  bed  at  bed  des- 

destruction 

truction 

Settling  type 

Collective  subsidence 

Free  settling 

Type  of  flow 

Type  A 

Type  B 

curve 

Existence  of  Ty|  (yield 

No  Ty]  (yield  stress) 

stress) 

Smooth  flow  curve 

Sharp  bend  in  flow  curve 

(power-law  type) 

•From  Ohtsubo  and  Muraoka  (1986). 


secretions  of  algae,  bacteria,  diatoms,  or  other  members  of  the 
benthic  and  epiflora  and  fauna  common  in  coastal  and  oceanic 
environments  (McCave,  1984).  Even  though  neither  the  phys¬ 
icochemical  nor  the  biochemical  bonding  forces  can  be  simply 
predicted  or  easily  measured  at  the  present  time,  advances  have 
been  made  in  cohesive  sediment  dynamics  in  the  last  30  years. 
Major  contributions  were  made  by  Einstein  and  Krone  (1962), 
Krone  (1962,  1963,  1986),  Partheniades  (1962,  1965,  1977, 
1986),  Migniot  (1968),  McCave  (1970,  1984),  Mehta  (1973, 
1986),  Owen  (1977),  and  Hunt  (1982,  1986). 

More  recently,  Mimura  (1989)  presented  a  comprehensive 
review  of  cohesive  sediment  dynamics  studies  in  Japan.  Signifi¬ 
cant  experiments  on  the  critical  shear  stress  and  rate  of  resuspen¬ 
sion  were  carried  out  by  Ohtsubo  and  Muraoka  (1986).  Forty 
types  of  sediment  were  tested  including  clays  (kaolinite,  mont- 
morillonite,  and  bentonite)  and  muds  taken  from  different  lakes 
and  estuaries.  To  examine  the  oiganic  effects,  starch  and  agar 
were  added,  respectively,  to  the  clay,  and  some  of  the  natural 
muds  were  tested  chemically.  The  most  important  result  showed 
that  cohesive  sediment  could  be  classified  into  two  groups  by  the 
difference  in  the  resuspension  behavior,  settling  characteristics, 
and  rheological  properties  over  a  wide  range  of  sediment  types 
and  water  contents.  Kaolinite,  kaolinite  containing  oiganic  matter, 
montmorillonite,  and  natural  muds  were  classified  into  one 
group,  namely  the  first  mud  group.  Bentonite  and  bentonite  con¬ 
taining  organic  matter  were  classified  as  the  second  mud  group. 
Some  important  characteristics  of  the  two  mud  groups  are  listed 
in  Table  20.2.  The  two  different  settling  types,  collective  subsi- 


20.  Riverine  and  Continental  Shelf 


195 


Figure  20.1.  Shear  stress-strain  rate  relationship  for  two 
mud  groups  (A)  Type  A  has  a  yield  stress;  (B)  Type  B  has 
no  well-defined  yield  strength.  (From  Ohtsubo  and 
Muraoka,  1986.) 


(A) 


(B) 


dence  and  free  settling,  indicate  the  difference  in  cohesive  forces 
between  the  two  mud  groups.  In  the  first  mud  group,  collective 
subsidence  settling  occurs  as  the  suspended  particles  aggregate 
and  form  a  particular  chain-like  structure.  The  different  types  of 
relationship  between  shear  stress  and  strain  rate  were  obtained 
as  shown  in  Figure  20. 1 .  Type  A  curve,  which  has  a  yield  stress 
Ty,  and  a  sharp  bend  indicative  of  material  having  a  finite  shear 
strength,  was  obtained  for  the  first  mud  group.  Type  B  curve, 
which  has  no  well-defined  yield  strength,  was  representative  of 
the  second  mud  group.  The  existence  of  the  yield  stress  Tyl  sug¬ 
gests  that  the  first  mud  group  has  strong  bonding  and/or  struc¬ 
ture  between  particles. 

The  sizes  of  hydrated  ions  are  expected  to  play  an  important 
role  in  the  aggregation  of  suspended  particles.  Some  of  the 
important  possible  mechanisms  for  clay-saltwater  interactions 
include  (1)  adsorbed  hydrated  cations  which  form  the  diffuse 
double  layer,  (2)  hydrogen  bonding,  (3)  hydration  of  exchange¬ 
able  cations,  (4)  charged  surface-dipole  attraction,  (5)  attraction 
by  van  der  Waals  forces,  (6)  electroosmosis,  (7)  chemicoosmo- 
sis,  (8)  thermal  diffusion,  (9)  external  body  forces,  and  (10) 
hydrodynamic  forces. 


Flow  Dynamics 

Water  Mass  Dynamics  on  the  Continental  Shelf 

The  transport  of  sediment  in  the  water  column  and  as  bed  load 
material  is  strongly  influenced  by  both  the  fluid  motion  as  well 
as  the  microstructure  in  suspension.  The  dynamics  and  micro¬ 
scale  characteristics  of  the  water  masses  also  significantly  affect 
the  developmental  history  of  fine-particle  microstructure.  On 
the  continental  shelf,  a  wide  range  of  spatial  and  temporal  scales 
of  water  motion  exists.  Since  long-term  changes  involve  deposi¬ 
tion  and  low  rate  transport  of  clay  materials,  water  motion  of  low 
intensity  is  important  over  geological  time  scales.  Wave-induced 
water  motion  and  currents,  tidal-generated  currents,  and  cur¬ 
rents  induced  by  fluid  density  differences  are  important  for  sedi¬ 
ment  transport  under  certain  environmental  conditions.  A  brief 
discussion  of  the  various  aspects  of  fluid  flow  important  in  sedi¬ 


ment  transport  studies  was  presented  by  Weggel  (1972)  and  a 
review  of  the  coastal  processes  was  given  by  Johnson  and  Eagle- 
son  (1966).  According  to  Johnson  and  Eagleson,  ocean  currents 
can  be  classified  into  several  groups;  (1)  currents  generated  by 
density  distributions  in  the  sea,  (2)  currents  resulting  from  wind 
stresses,  (3)  tidal  currents  and  those  associated  with  internal 
waves,  (4)  currents  caused  by  the  activity  of  surface  waves,  and 
(5)  currents  induced  by  the  intrusion  of  freshwater  entering  the 
ocean  at  river  mouths.  The  various  types  of  currents  are  depicted 
in  Figure  20.2. 

The  flow  regions  in  a  vertical  section  of  continental  shelf 
water  are  shown  schematically  in  Figure  20.3.  The  flow  charac¬ 
teristics  for  the  “wall”  region  are  very  significant  in  sediment 
transport  and  the  dynamics  are  different  from  the  main  water 
mass  body.  The  turbulence  and  dynamics  of  the  “wall”  region  are 
of  most  importance  in  terms  of  sediment  transport  and  sus¬ 
pended  sediment  microstructure. 


Turbulence 

Turbulence  is  the  irregular  motion  of  rotational  flows.  The  ther¬ 
modynamic  and  hydrodynamic  variables  of  these  flows  depict 
the  chaotic  fluctuations  in  space  and  time  coordinates.  The  fol¬ 
lowing  possible  mechanisms  of  turbulence  generation  in  the 
ocean  were  pointed  out  by  Monnin  and  Ozmidow  (1985); 

1.  Instability  of  vertical  velocity  gradients  in  drift  flow-these 
velocity  gradients  are  induced  by  direct  effects  of  the  wind  on 
the  upper  mixed  layer  (UML). 

2.  Overturning  of  surface  waves -this  mechanism,  and  to  a  less 
extent  the  hydrodynamic  instability  of  wave  motions  induced 
by  surface  waves,  is  the  most  powerful  source  of  UML  turbu¬ 
lence. 

3.  Instability  of  vertical  velocity  gradients  in  stratified  large- 
scale  oceanic  flows -because  of  the  large  scale  and  low  rate  of 
turbulent  energy  generation,  this  mechanism  is  considered 
not  important  for  shelf  sediment  transport. 

4.  Hydrodynamic  instability  of  quasihorizontal  mesoscale  non¬ 
stationary  flows -empirical  evidence  exists  for  the  genera¬ 
tion  of  small-scale  turbulence. 
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Figure  20.2.  Classification  of  ocean  currents.  (From  John¬ 
son  and  Eagleson,  1966.) 


5.  Instability  of  local  velocity  gradients  in  internal  waves-the 
instability  of  internal  waves  in  the  turbulent  microstructure  of 
the  water  increases  the  turbulent  energy  by  a  very  small 
amount,  which  results  in  only  slight  dissipation. 

6.  Convection  in  layers  with  unstable  density  stratification - 
this  mechanism  can  create  turbulent  (microstructure)  layers. 


^Sea  surface 

upper  mixed  layer 


Bulk  ocean 


•Bottom  boundary  layer 
'Wall'  region  _ 


Figure  20.3.  Flow  regions  in  a  vertical  section  (not  to  scale). 


7.  Instability  of  vertical  velocity  gradients  in  the  bottom  bound¬ 
ary  layer  (BBL)-this  type  of  turbulence  is  created  by  Ekman 
layer  instability. 

In  addition  to  the  turbulence  generation  mechanisms  listed 
above,  the  dynamics  and  turbulence  characteristics  in  the  “wall” 
region  are  important  considerations  in  continental  shelf  water 
dynamics.  The  thickness  of  the  “wall”  region  (Nychas  et  al., 
1973)  is  defined  by 

8W  =  70  v/h.  (1) 

where  u,  (u.  =  tb1/2p'l/2,  here  xh  is  the  bed  shear  stress  and  p  is 
the  density  of  the  water)  is  the  friction  velocity  at  the  wall,  and  v 
is  the  kinematic  viscosity  of  the  water.  Dirbulence  in  the  “wall” 
region  is  not  isotropic  and  not  homogeneous.  The  turbulence, 
however,  is  termed  isotropic  if  its  statistical  properties  have  no 
preferential  direction,  and  it  is  said  to  be  homogeneous  if  a  ran¬ 
dom  motion  exists  whose  average  properties  are  independent  of 
position  in  the  fluid.  Classical  works  of  turbulence  were  pre¬ 
sented  by  Hinze  (1959),  Batchelor  (1960),  and  Townsend  (1976). 


Dynamics  in  the  “Wall”  Region 

A  thin  viscous  sublayer  exists  next  to  the  interface  in  the  “wall” 
region  for  a  smooth,  rigid  boundary.  The  ranges  of  the  viscous 
sublayer  and  the  “wall"  region  can  be  estimated  as  follows: 
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Table  20.3.  Estimated  thickness  ranges  for  viscous  sub¬ 
layer  and  “wall"  region. 

Friction  velocity  u„  (cm/sec) 

0.1  1.0  10 

Viscous  sublayer  (cm)  <0.5  <0.05  <0.005 

“Wall”  region  (cm)  <7.0  <0.7  <0.07 

Viscous  sublayer  y  <  5 v/u.  (2) 

“Wall"  region  y  <  70 v/«.  (3) 

where  u.  is  the  friction  velocity  and  v  is  the  kinematic  viscosity 
of  the  water.  Assuming  v  =  0.01  cm2/sec,  the  thickness  ranges 
are  estimated  in  Table  20.3.  The  condition  for  a  boundary  to  be 
hydraulically  smooth  is  that  the  dimensionless  roughness  Rey¬ 
nolds  number 

u.k/v  <  5  (4) 

where  k  is  the  roughness  parameter  (Schlichting,  1960).  Accord¬ 
ing  to  Einstein  ( 1950).  the  total  friction  velocity  (or  shear  stress) 
can  be  expressed  as 

=  u,  +  u,  (  (5) 

where  u.  is  the  "pure  friction”  component  interacting  between 
the  fluid  and  the  bed  surface,  irrespective  of  whether  the  surface 
is  a  flat  plane  or  irregular,  and  u.,  is  the  resultant  of  the  drag 
forces  acting  on  the  irregularities.  Discussions  of  the  bottom 
friction  were  given  by  Yalin  (1977).  Based  on  Einstein’s  sugges¬ 
tion  as  stated  in  Eq.  (5),  the  sea  bottom  is  likely  to  be  classified 
as  smooth.  Even  for  the  rough  boundary  (5  <  u.k/v  <  70, 
where  k  is  the  grain  roughness  parameter),  the  concept  of  the 
“wall”  region  is  still  valid. 

The  high  shear  and  the  dynamics  of  the  viscous  sublayer  and 
the  “wall”  region  are  the  key  mechanisms  that  play  a  role  in  the 
deposition,  resuspension,  and  floe  breakup.  From  turbulent  struc¬ 
ture  measurements  in  flat  plate  and  pipe  flow,  the  regions  of  maxi¬ 
mum  production  and  maximum  dissipation  are  just  outside  the 
viscous  sublayer  (Hinze,  1959;  Townsend,  1976).  It  is  important 
to  note  that  a  dynamic  viscous  sublayer  model  was  proposed  by 
Einstein  and  Li  (1955,  1956,  1958).  The  proposed  model  predicts 
a  periodic  growth  and  decay  of  the  viscous  sublayer.  As  a  conse¬ 
quence  of  this  model  some  estimates  can  be  made  for  the  genera¬ 
tion  of  turbulence  near  the  boundary.  The  unsteady  and  intermit¬ 
tent  model  was  based  on  visual  observations,  pressure  fluctuation 
measurements,  and  logic.  Observations  revealed  that  dye  and  sed¬ 
iment  particles  were  exchanged  between  a  smooth  boundary  and 
the  turbulent  flow  through  the  viscous  sublayer,  thus  making  the 
assumption  of  a  quasisteady  sublayer  impossible  (Einstein  and  Li, 
1958).  Additional  visual  and  quantitative  studies  of  the  viscous 
sublayer  and  “wall”  region  were  carried  out  by  Kline  et  al.  (1967). 
An  experimental  investigation  with  hot-wire  anemometers  was 
performed  by  Gupta  et  al.  (1971)  to  evaluate  the  spatial  structure 


in  and  near  the  viscous  sublayer.  Later  experiments  (Nychas  et 
al.,  1973;  Offen  and  Kline,  1975)  confirmed  that  the  production 
of  the  Reynolds  stress  is  intermittent  and  is  associated  with  a 
sequence  of  motions  called  “bursting”  (Offen  and  Kline,  1974, 
1975;  Dyer,  1986). 


Wave  Boundary  Layer 

In  continental  shelf  waters,  the  wave  action  may  reach  the  sea 
bottom.  It  is  well  known  that  the  irrotational  theory  can  be 
applied  to  the  entire  gravitational  wave  motion  except  for  a  very 
thin  layer  adjacent  to  solid  boundaries.  Near  solid  boundaries 
the  viscous  effect  cannot  be  neglected  as  compared  with  the  iner¬ 
tia  forces,  and  the  thin  layer  can  be  treated  according  to  the 
boundary  layer  theory.  An  experimental  study  was  carried  out 
by  Li  (1954)  to  investigate  the  wave  boundary  stability  by  oscil¬ 
lating  a  plate  in  still  water  with  the  range,  amplitude,  and  period 
as  computed  by  wave  theory.  According  to  Li  (1954),  the  thick¬ 
ness  of  the  occillatory  boundary  layer  is  defined  as 

8,  =  6.5  y/vTw  (6) 

where  v  is  the  kinematic  viscosity  and  to  is  the  angular  fre¬ 
quency.  The  thickness  of  the  occillatory  boundary  layer  is  very 
thin.  For  an  occillatory  period  of  10  sec,  8,  has  the  thickness 
about  0.82  cm. 

The  critical  Reynolds  number  at  which  the  transition  from 
laminar  to  turbulent  takes  place  is  found  to  be  a  constant 

R  =  G)l,2d,/vl/2  =  800  (7) 

where  d,  is  the  maximum  amplitude  of  horizontal  displacement 
at  the  boundary.  The  critical  Reynolds  number  should  be  lower 
in  the  field  (in  situ )  than  the  laboratory  values  stated  above  since 
the  flow  disturbances  are  much  higher  in  the  field. 

Further  information  for  wave  boundary  layers  and  their  rela¬ 
tion  to  sediment  transport  was  provided  by  Teleki  (1972)  and 
Einstein  (1972).  Laboratory  experiments  were  carried  out  by 
Kemp  and  Simons  (1982)  to  investigate  the  interaction  between 
waves  and  currents.  Theoretical  models  were  presented  by 
Smith  (1977),  Grant  and  Madsen  (1979),  Grant  et  al.  (1984), 
Herbich  et  al.  (1984),  and  Shen  (1986). 


Detrital  Particles  and  Flow  Interactions 
Fluid-Particle  Dynamics 

The  transport  of  solid  particles  by  random  fluid  motions  that 
may  include  molecular  motion  and  turbulence  is  diffusive  in 
nature.  For  cohesive  particles,  the  particle-particle  interaction 
and  the  change  of  particle  sizes  make  the  problem  very  complex. 
Several  scales  of  interaction  are  involved  in  the  particle  and  flow 
dynamics;  molecular  scale,  solid  particle  size,  turbulent  eddy 
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Ibble  20.4.  Densities  and  shear  strengths  of  aggregates  in  the  capillary 
viscometer* 


Sample 

Density 

(g/cm3) 

Shear  strength 
(dyn/cm2) 

Cation  exchangeable  capacity3 
(meq/100  g) 

Wilmington  District 

1.250 

21.0 

32 

Brunswick  Harbor 

1.164 

33.7 

38 

Gulfport  Channel 

1.205 

45.9 

49 

San  Francisco  Bay 

1.269 

21.8 

34 

Delaware  River 

1.305 

11 

23 

Potomac  River 

1.437 

19 

15 

White  River 

1.212 

48.6 

60 

*From  Krone  (1986). 

3  Cation  exchange  capacity  of  the  <  10  pm  fraction. 


size,  and  magnitude  of  the  mean  motion.  Because  of  the  random 
nature  of  sediment-fluid  dynamics  and  the  lack  of  complete 
understanding  of  the  problems  and  microscale  events,  both 
statistical  and  semiempirical  approaches  are  required  to 
approach  a  solution  to  the  sediment  transport  processes. 

For  clay  particles  to  aggregate,  mechanisms  must  exist  to 
bring  particles  into  collisions  and  the  collisions  must  result  in 
flocculation  (coagulation).  The  flocculation  process  is  due  to  the 
physicochemical  interaction  of  the  clay  surfaces  and  surround¬ 
ing  electrolytic  fluid.  Common  physical  mechanisms  involving 
particle  collisions  are  Brownian  motion,  laminar  and  turbulent 
shear,  and  differential  settling  (see  Bennett  et  al.,  this  volume). 
A  general  equation  for  the  particles  undergoing  flocculation 
including  their  size  distribution  was  given  by  Friedlander 
(1977).  The  “coagulation  kernels”  were  found  for  the  Brownian 
motion,  laminar  shear,  turbulent  shear,  and  differential  settling. 
Analogous  to  the  work  in  the  universal  equilibrium  concept  of 
turbulence,  a  similar  solution  for  coagulation  was  developed  by 
Friedlander  (1960)  and  later  extended  by  Hunt  (1982).  Similar 
to  flocculation,  the  equation  for  the  particle  size  distribution 
undergoing  breakup  was  given  by  Jeffrey  (1982);  however,  the 
“breakup  kernels”  were  not  found.  Additional  data  on  the  coagu¬ 
lation  and  breakup  were  presented  by  McCave  (1984)  and  Hunt 
(1986).  Based  on  known  measurements,  the  size  distribution  can 
be  estimated  for  shelf  environments  when  the  system  is  free  of 
biological  influences:  (1)  particle  sized  spectra  are  flat  at  low 
concentrations  (<  1  x  10-5  g/cm3)  in  the  size  range  1-60  pm 
(Eisma  and  Kalf,  1979),  and  (2)  at  moderate  concentration  (~ 

1  x  10 _4  g/cm3),  a  modal  peak  occurs  between  5  and  20  pm 
(Kranck,  1973,  1981).  Several  modal  peaks,  however,  may  exist 
in  the  shelf  where  plankton  are  important.  Once  the  size  of  the 
floe  is  known,  the  settling  velocity  of  a  single  floe  can  be  esti¬ 
mated  theoretically  or  empirically. 

Krone  (1986)  presented  an  hypothesis  that  described  the 
orders  of  aggregates  and  laboratory  measurement  results  of  the 
capillary  viscometer  for  density,  shear  strengths  of  aggregates, 
and  cation  exchange  capacity  (Table  20.4).  The  Bingham  shear 
strength  correlates  with  the  cation  exchange  capacity  as 

(8) 


Table  20.5.  Properties  of  sediment  aggregates* 


Sediment  sample 

Order  of 
aggregation 

Density3 

(g/cm3) 

Shear  strength 
(dyn/cm2) 

Porosity5 

(«) 

Wilmington  District 

0 

1.250 

21 

0.862 

1 

1.132 

9.4 

0.934 

2 

1.093 

2.6 

0.958 

3 

1.074 

1.2 

0.970 

Brunswick  Harbor 

0 

1.164 

34 

0.914 

1 

1.090 

4.1 

0.960 

2 

1.067 

1.2 

0.974 

3 

1.056 

0.62 

0.981 

Gulfport  Channel 

0 

1.205 

46 

0.889 

1 

1.106 

6.9 

0.950 

2 

1.078 

4.7 

0.967 

3 

1.065 

1.8 

0.975 

San  Francisco  Bay 

0 

1.269 

22 

0.850 

1 

1.179 

3.9 

0.905 

2 

1.137 

1.4 

0.931 

3 

1.113 

1.4 

0.966 

4 

1.098 

0.82 

0.955 

5 

1.087 

0.36 

0.962 

6 

1.079 

0.20 

0.967 

White  River 

0 

1.212 

49 

0.885 

(in  salt  water) 

1 

1.109 

6.8 

0.948 

2 

1.079 

4.7 

0.967 

3 

1.065 

1.9 

0.975 

♦From  Krone  (1986). 

3  Density  in  sea  water  =  1.025  g/cm’. 
3  Solid  particle  density  =  2.65  g/cm5. 


In  Table  20.5,  the  results  of  the  rotating  cylinder  viscometer  are 
shown  for  the  order  of  aggregation,  density,  and  shear  strength, 
as  well  as  the  porosity.  Krone  stated  that 

1 .  The  order  of  the  suspended  aggregates  depends  on  the  history 
of  the  suspension:  such  as  if  the  particles  have  been  sus¬ 
pended  in  fluid  zones  having  very  high  velocity  gradients  or 
if  they  are  newly  formed  from  dispersed  mineral  particles, 
the  order  would  be  zero.  Aggregates  that  have  formed  in 
moderate  velocity  gradients  were  classified  as  one  order 
higher. 

2.  Settling  aggregates  that  adhere  to  the  bed  surface  during 
deposition  have  a  structure  that  at  the  bed  surface  is  classified 
as  one  order  higher  than  that  of  the  depositing  aggregate.  The 
critical  shear  stress  for  deposition  (rcd)  also  is  considered 
one  order  of  aggregation  higher  than  the  order  of  the  deposit¬ 
ing  aggregate. 

Flume  studies  of  deposition  of  suspended  cohesive  sediment 
(Krone,  1986)  showed  that 

dcldt  =  ~{Wsdh)  (1  -  T„/Tcd)  (9) 

where  c  is  the  concentration,  t  is  the  time,  1TS  is  the  settling 
velocity,  h  is  the  water  depth,  rh  is  the  bed  shear  stress,  and  rtd 
is  the  critical  shear  stress  for  deposition.  The  term  (1  -  Tb/rcd) 
expresses  the  probability  of  the  settling  aggregates  to  adhere  to 
the  bed. 


tb  =  3.92  +  0.8447  (CEC) 
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Table  20.6.  Values  of  rcd 

for  three  sediments.* 

Sediment 

•tcd  (dyn/cm2) 

Kaolinite 

1.5 

Bay  mud 

1.0 

Maracaibo  mud 

0.8 

*From  Mehta  (1986). 


Mehta  (1973,  1986)  carried  out  extensive  deposition  tests 
under  steady  flows  and  reanalyzed  some  previous  flume  data. 
He  found  that  the  critical  shear  stress  xcd  depended  on  the  sedi¬ 
ment  composition,  and  varied  from  0.4  to  1.5  dyn/cm2  The 
values  of  xcd  are  shown  in  Table  20.6  for  three  sediments. 

Partheniades  (1962,  1965,  1986)  conducted  experimental  and 
theoretical  studies  for  both  erosion  and  deposition.  The  first 
experiments  were  conducted  on  a  remolded  bed  at  a  comparable 
field  density  and  of  uniform  consistency  about  1  in.  thick.  At 
constant  bed  shear  stresses  above  a  threshold  value,  rce  =  0.96 
dyn/cm2,  the  rates  of  erosion  were  constant  and  independent  of 
the  suspended  sediment  concentration  (Partheniades,  1962, 
1965,  1986).  Data  and  observations  indicated  that  no  simultane¬ 
ous  erosion  and  deposition  of  sediment  occurred  during  the 
experiments.  Based  on  statistical  considerations,  models  were 
developed  by  Partheniades  for  erosion,  deposition,  and  sedi¬ 
ment-flow  interactions.  The  stochastic  model  for  sediment-flow 
interaction  was  based  on  an  extension  of  the  Einstein  (1950) 
approach  for  the  bed-load  function.  The  model  indicates  that 
simultaneous  erosion  and  deposition  of  sediment  particles  may 
occur  only  in  a  limited  range  for  cohesive  sediments.  Laboratory 
studies  of  erosion  and  deposition  are  very  useful  for  studying 
and  describing  shelf  sediment  transport  processes  and  clay 
microstructure  development. 


Turbulent  Scale ,  Dissipation  Rate,  and  Mean  Shear  Rate 

Let  us  now  address  the  subject  concerning  the  scale  of  turbulence. 
For  isotropic  turbulence,  Kolmogoroff  made  the  following 
hypothesis  (Hinze,  1959):  “At  sufficiently  high  Reynolds  numbers 
there  is  a  range  of  high  wave-numbers  where  the  turbulence  is 
statistically  in  equilibrium  and  uniquely  determined  by  the 
parameters  c  and  v.  This  state  of  equilibrium  is  universal”  (Hinze, 
1959,  p.  184).  e  and  v  describe  the  rate  of  turbulent  energy  dissi¬ 
pation  and  the  molecular  kinematic  viscosity,  respectively. 

From  dimensional  reasoning,  the  characteristic  scales  can  be 
obtained  from  the  equilibrium  range  as 

Length  scale:  r|  =  (v7e)i/4  (10) 

Velocity  scale:  u  =  (ve)l/4  (11) 

Shear  rate  scale:  G  =  (e/v),/2  (12) 

The  rate  of  turbulent  energy  dissipation  (or  production),  e,  in  the 
ocean  was  estimated  by  Monin  and  Ozmidov  (1985)  as  (1)  in 


Ibbte  20.7.  Characteristic  scales  in  equilibrium  range  for  three  layers  of  the 
ocean. 


Region 

Dissipation 

rate 

e  =  (cm2/secJ) 

Length 

scale 

0  (cm) 

Velocity 

scale 

o  (cm/sec) 

Shear  rate 
scale 

C  (sec"') 

Upper  mixed 

I  O'2- 10-' 

I0>-0.56 

1  O’ 1 - 1 . 8 

1-3.2 

layer 

x  I0'1 

x  10-' 

Bulk  ocean 

io-5 

5.6  X 

1.8  x 

3.2  x 

io-' 

IO’2 

IO'2 

Bottom  boundary 

to-* 

3.2  x 

3.2  x 

IO"' 

layer 

10-' 

IO'2 

upper  mixed  layer,  8  =  10'2  -  10"1  cm2/sec3;  (2)  in  bulk  ocean, 
e  =  10's  cm2/sec3;  and  (3)  in  BBL,  with  turbulence  created 
chiefly  by  Ekman  layer  instability,  £  =  10*4  cm2/sec3.  Based  on 
these  estimates,  the  characteristic  scales  were  determined  (Table 
20.7).  It  is  important  to  note  that  the  assumption  of  local 
isotropy  and  equilibrium  range  is  reasonable  in  all  the  three 
layers  except  in  the  “wall”  region.  In  the  “wall”  region,  the  rate 
of  energy  dissipation  and  the  mean  shear  rate  are  estimated 
based  on  the  following  relationships: 

In  the  viscous  sublayer,  the  viscous  dissipation  rate  is 

£v  =  m.4/v  (13) 

and  mean  shear  rate  is 

du/dy  =  m.2/v  (14) 

At  the  edge  of  the  “wall”  region  the  turbulent  dissipation  rate  is 
E  =  M,3/Ky  (15) 

and  mean  shear  rate  is 

du/dy  =  K./icy  (16) 

where  k  is  the  von  Karman’s  constant  and  has  the  value  of  0.4  for 
clean  fluids  (may  be  lower  with  sediment). 

The  estimates  of  flow  characteristics  are  given  in  Tables  20.8 
and  20.9  for  the  viscous  sublayer  and  for  the  edge  of  the  “wall" 
region.  When  the  friction  velocity  reaches  the  value  of  1.0 
cm/sec,  high  viscosity  dissipation  and  high  mean  shear  rates 
occur  in  the  viscous  sublayer  (Table  20.8)  and  the  turbulent  dis¬ 
sipation  rates  at  the  edge  of  the  “wall”  region  (Table  20.9)  are 
higher  than  that  of  the  dissipation  rate  in  the  upper  mixed  layer. 


Table  20.8.  Eslimale  of  flow  characteristics  in  the  viscous  sublayer. 


Friction  velocity  u,  (cm/sec) 


0.1 

1.0 

10 

Viscous  dissipation  rate 

IO'2 

102 

IO* 

ev  (cm2/sec3) 

Mean  shear  rate 

1 

I02 

10* 

du 

.  sec  ') 
dy 
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Table  20.9.  Estimate  of  flow  characteristics  at  the  edge  of  “wall"  region. 


Friction  velocity  u,  (cm/sec) 


At  v  =  70  u/w. 

0.1 

1.0 

10 

Turbulent  dissipation  rate 

e  (cm2/sec’t 

3.6  x  10-“ 

3.6 

3.6  x  104 

Mean  shear  rate 

u  (sec1) 
dy 

3.6  x  10"2 

3.6 

3.6  x  I02 

. 

_ 

_ _ 

bulk  ocean,  and  bottom  boundary  layer.  Due  to  the  fact  that  the 
turbulent  dissipation  rate  and  the  turbulent  production  rate  are 
much  higher  inside  the  “wall”  region  than  outside,  and  because 
of  the  unsteady,  intermittent,  dynamic  nature  of  this  region,  the 
breakup  of  floes  takes  place  mainly  in  the  "wait'  region.  On  the 
other  hand,  the  floes  also  may  affect  the  turbulent  characteris¬ 
tics,  and  the  structure  of  the  “wall”  region,  the  fluid  and  solids 
interacting  synergistically.  The  turbulence  in  the  other  layers 
contributes  to  the  mixing,  flocculation,  diffusion,  and  deforma¬ 
tion  of  the  floes,  but  at  generally  lower  dynamic  energies  than  in 
the  “wall”  region. 

Description  of  the  Particles  and  Flow  Interactions  on  the  Shelf 
and  Riverine  Environments:  Significance  to  Clay  Microstructure 

Despite  the  complexity  of  the  particle  types  and  morphology  and 
flow  dynamics  on  the  continental  shelf,  the  mechanics  of  solid 
particle  or  floe  movement  can  be  described  and  approximated  as 
follows: 

1 .  The  motion  of  any  solid  particle  (or  aggregate)  is  determined 
by  Newton’s  laws.  The  forces  acting  on  the  particle  include 
hydrodynamic  forces,  electrochemical  particle-particle  inter¬ 
action  forces,  and  gravitational  force.  Due  to  the  various 
shapes  of  particles  and  the  random  nature  of  the  force 
actions,  the  particle  may  reorientate  and  floccules  may 
deform. 

2.  As  an  example,  suppose  we  follow  a  solid  particle  (aggregate) 
in  a  location  outside  the  “wall”  region  (i.e.,  yu./v  >  70). 
Before  a  particle-particle  collision,  the  motion  of  this  parti¬ 
cle  can  be  considered  as  a  single  particle  that  is  subjected  to 
the  actions  of  the  steady  current,  periodic  motions  such  as 
wave  and  tide,  flow  turbulence.  Brownian  motion,  and 
gravitational  force.  This  particle  reorientates  or  moves  more 

or  less  in  a  random  manner.  However,  the  motion  of  the  parti¬ 
cle  has  the  tendency  to  follow  a  downward  trajectory  due  to 
the  gravitational  force.  Once  particle-to-particle  collision 

occurs,  the  compound  particle  may  either  bond  or  separate 
depending  on  the  impact  force  and  the  physicochemical 
forces.  The  strength  of  the  newly  formed  aggregate  is  proba¬ 
bly  dependent  on  (a)  the  number  of  molecular  bonds  at  the 


particle-particle  contact,  (b)  the  type  of  available  hydrated 
ions  in  the  double  layer,  (c)  the  presence  or  absence  of  other 
molecules  (organic  and  inorganic)  in  proximity  to  the  points 
of  contact,  and  (d)  the  geometry  of  the  contact  area  and  rela¬ 
tive  orientation  of  the  particles  [edge-to-edge  (E-E)  and  edge- 
to-face  (E-F)  contact  generally  would  be  weaker  than  face-to- 
face  (F-F)  contacts].  After  the  particle-particle  collision,  the 
new  aggregate  will  be  considered  as  a  single  particle  as  before 
and  the  processes  will  continue  in  the  flow  regime  outside  the 
"wall”  region. 

3.  As  the  aggregate  drops  inside  the  “wall”  region  and  into  the 
viscous  sublayer,  the  particle  may  deform,  break  up,  and  eject 
outside  the  region  due  to  the  high  shear  and  dynamics  of  this 
region.  Because  of  the  unsteady  and  intermittent  nature  in  the 
“wall”  region,  the  high  lift  force  may  prevent  aggregate  depo¬ 
sition.  However,  the  aggregate  may  pass  through  the  “wall” 
region  and  be  deposited  at  the  interface  if  the  bottom  friction 
velocity  is  below  a  critical  value  (mv, i  —  0.6-1. 3  cm/sec). 
Simultaneous  erosion  and  deposition  generally  do  not  occur. 
The  deposited  sediment  will  consolidate,  be  reworked,  and  be 
altered  by  the  mechanical,  chemical,  and  biological  forces 
and  mechanisms  driving  microfabric  development  (Bennett 
et  al.,  this  volume).  As  the  bottom  friction  velocity  reaches 
another  critical  value  («^.e  =  0.6-3. 0  cm/sec)  at  another 
time,  the  deposited  sediment  will  be  resuspended. 

4.  Since  the  nature  of  the  problem  has  virtually  infinite  degrees 
of  freedom,  the  high  variability  in  time,  space,  and  scale,  both 
a  statistical  concept  and  semiempirical  approaches  are 
required  to  develop  meaningful  predictive  models  of  sedi¬ 
ment  transport  and  deposition  and  the  processes  and  mechan¬ 
isms  in  microfabric  development. 

Based  on  this  study  and  the  above  descriptions  and  the  discus¬ 
sions  in  previous  sections,  the  suspended  clay  microstructure, 
can  be  described  as  follows:  (1)  clay  mineral  particle  microstruc¬ 
ture  (fabric  and  physicochemistry)  are  the  fundamental  “build¬ 
ing  blocks”;  (2)  the  available  particle  size  distribution  affects  the 
aggregate  size  distribution,  which  commonly  is  flat  in  the  range 
1-50  pm  (where  plankton  are  important  several  modes  can 
exist);  (3)  the  density  of  the  aggregates  have  the  range  1 .06-1 .4 
g/cm3;  (4)  after  the  initial  contact,  which  is  dominantly  edge-to- 
face  in  the  marine  environment,  stairstep  structure  (stronger  and 
more  stable  structure)  is  preferred  due  to  the  high  shear  in  the 
viscous  sublayer  and  the  “wall"  region;  and  (5)  chain-like  and 
vortex-like  structures  may  exist  partially  due  to  turbulence  espe¬ 
cially  in  the  “wall”  region. 


Summary  and  Conclusions 

The  suspended  clay  particles  in  salt  water  may  coagulate  due  to 
the  effects  of  electrostatic  forces,  van  der  Waals  forces,  and  bio¬ 
chemical  activities.  The  flocculation-deflocculation  behavior  of 
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clays  in  suspension  is  dependent  on  many  physicochemical 
properties,  especially  the  clay  mineral  types,  their  particle  size 
distribution,  and  the  ionic  characteristics  of  the  fluid.  The  parti¬ 
cles  and  flow  interactions  such  as  Brownian  motion,  laminar  and 
turbulent  shear,  and  differential  settling  are  of  importance  in  the 
flocculation  behavior. 

A  thin  water  layer  near  the  bottom  plays  a  very  important  role 
on  the  suspended  clay  particle  and  aggregates  and  on  the  parti¬ 
cles  at  the  depositional  interface.  This  thin  layer  is  called  a 
“wall”  region  and  has  a  thickness  equal  to  70v/m,,  where  v  is  the 
kinematic  viscosity  of  the  fluid  and  u ,  is  the  “pure”  friction 
velocity.  Outside  the  “wall”  region,  the  turbulence  and  the  mean 
shear  play  a  role  in  flocculation.  Inside  the  “wall”  region,  the 
particle  aggregate  may  break  up  and  be  ejected  outside  the 
region  due  to  the  high  shear  and  dynamics  of  the  “wall”  region 
and  the  microfabric  of  the  floccules  may  be  changed.  The 
aggregate  may  pass  through  the  “wall”  region  and  be  deposited 
on  the  bottom  if  the  friction  velocity  is  below  a  critical  value 
(e.g.,  wVd  =  0.6-1 .3  cm/sec). 

The  friction  velocity  due  to  the  drag  forces  acting  on  the 
irregularities,  u.,,  is  a  contribution  mainly  outside  the  “wall” 
region  insofar  as  the  particle-flow  interaction  is  concerned.  Due 
to  the  distinct  roles  of  the  u.  and  u,„  the  application  of  the 
logarithmic  law  is  not  sufficient  for  the  determination  of  the 
friction  velocity  and  the  roughness  length. 

The  interaction  of  waves  and  currents  is  complex.  For  the  fric¬ 
tion  velocity,  a  linear  superposition  is  recommended  as  the  first 
approximation. 
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